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Abstract 

We use the replica method in order to obtain an expression for the 
variational free energy of an Ising ferromagnet on a Viana-Bray lattice 
in the presence of random external fields. Introducing a global order 
parameter, in the replica-symmetric context, the problem is reduced 
to the analysis of the solutions of a nonlinear integral equation. At 
zero temperature, and under some restrictions on the form of the 
random fields, we are able to perform a detailed analysis of stability 
of the replica-symmetric solutions. In contrast to the behaviour of 
the Sherrington-Kirkpatrick model for a spin glass in a uniform field, 
the paramagnetic solution is fully stable in a sufficiently large random 
field. 
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1 Introduction 



Magnetic systems with quenched disorder, including spin glasses and fer- 
romagnets in a random field, have been intensively studied during the last 
decades pQ. There are many applications of these disordered systems, rang- 
ing from the study of the behavior of random magnets, which is a traditional 
ground test for the ideas of statistical mechanics, to the analysis of different 
sorts of optimization problems in distinct areas of science. The mean-field 
version of an Ising spin glass, with Gaussian distribution of exchange interac- 
tions, also known as the Sherrington-Kirkpatrick model, which can be solved 
by the replica method, displays a low-temperature glassy phase, character- 
ized by the instability of a replica-symmetric solution, which indicates the 
need of breaking replica symmetry and the existence of many ultrametrically 
organized states. In contrast to this rich behavior, the mean-field, Curie- 
Weiss, version of an Ising ferromagnet in a random field (RFIM), which can 
be solved without recourse to the replica method, leads to rather uninter- 
esting, replica-symmetric, exact solutions. There have been, however, some 
indications that a ferromagnet in a random field may have a glassy behavior. 

We were motivated to look again at this problem by a number of early 
and some more recent investigations of the RFIM. De Almeida and Bruinsma 
[2] have done some calculations, beyond the usual Curie- Weiss, mean-field, 
approximation, for analyzing the behavior of a bond-diluted Ising antiferro- 
magnet in a field. For large dimensionality, these calculations lead to the 
presence of a glassy region in the applied field versus temperature phase di- 
agram, between paramagnetic and ordered phases, which can be shown to 
disappear in the limit of infinite dimension. On the basis of the equivalence, 
at a mean-field level, between the critical behavior of a ferromagnet in a 
random field and of a dilute antiferromagnet in a uniform field, this result 
gives an indication of the possible existence of a glassy phase in the RFIM. 
A glassy behavior is also present in a recent "extended mean-field" calcu- 
lation by Pastor and collaborators 3j for the phase diagram of the RFIM. 
These results are claimed to agree with earlier work of Mezard and Young 
13] using a screening approximation in order to characterize the instability of 
the replica-symmetric solutions, to lowest order in 1/m, in a renormalization- 
group calculation for an m-component spin ferromagnet in a random field. It 
should mentioned that calculations for the RFIM on a Bethe lattice already 
indicate a rich ground-state structure [S] and peculiar hysteresis effects 0. 
Also, field-theoretical renormalizat ion-group calculations for a soft-spin ver- 
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sion of the RFIM, which were confirmed by a formulation of the dynamics, 
have shown the need to include extra terms involving replicas with different 
indices, which in turn may lead to an instability of the replica-symmetric 
solution in the paramagnetic region [7j. The prediction of a glassy phase in 
the calculations beyond the usual mean-field approximation were the main 
motivation to revisit this problem. We then decided to use a model devised 
by Viana and Bray jH], which is designed to gauge the effects of the (finite) 
connectivity of a lattice. 

According to the original work of Viana and Bray, we consider an Ising 
model with pair interactions Jjj, between all sites i and j, such that Jij = 
J > 0, with probability c/N, where N is the total number of sites, and 
Jij = 0, with probability 1 — c/N. This choice of interactions gives rise 
to the so-called Viana-Bray lattice. The parameter c > can be regarded 
as the (finite) mean connectivity per site. The solutions of an Ising spin- 
glass on the Viana-Bray lattice have been first analyzed in the vicinity of 
the transition temperature |H] . The more involved low-temperature behavior 
has been considered by Kanter and Sompolinsky [5], and Mezard and Parisi 
[TUj . The analysis of stability of the replica-symmetric solutions near perco- 
lation (c = 1) has been carried out by de Dominicis and collaborators [TT] . 
Although there are calculations for the RFIM on a Bethe lattice, there is no 
comprehensive analysis on a graph as the Viana-Bray lattice. According to 
the previous work for the Ising spin glass, we introduce a global order param- 
eter and formulate the replica-symmetric solutions for this problem in terms 
of an integral equation. In the ground state, and under some conditions on 
the form of the random fields, we show that the replica-symmetric solutions 
can be written as a series expansion. We have been able to perform a de- 
tailed analysis of stability of this replica-symmetric solution. In a sufficiently 
strong random field, the explicit calculation of the eigenvalues of a functional 
Hessian form shows the stability of the paramagnetic solution, which seems 
to preclude the existence of a glassy phase (in contrast to earlier expectations 
and the results of Pastor and collaborators 3j). Some numerical calculations 
confirm these findings for the paramagnetic solution, and indicate that the 
replica-symmetric ferromagnetic solution is also stable in the ferromagnetic 
region of the phase diagram. 

The layout of this paper is as follows. In Section 2, we define the model 
and formulate the replica-symmetric solutions in terms of an integral equation 
for the global order parameter. The analysis of stability of the paramagnetic 
solution is reported in Section 3. Some conclusions, and connections with 
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recent work, are presented in Section 4. 



2 Formulation of the problem 



The ferromagnetic Ising model on the Viana-Bray lattice is given by the 
Hamiltonian 



N 



i<Ji 



(1) 



(ij) i=l 

where {ij) refers to a pair of sites, o"i = ±1 for all sites, and {Jij} and {Hi} 
are sets of independent, identically distributed, random variables, associated 
with probability distributions 



and 



pj (Jij) = ^5 {Jij - J) + (^1 - ^) 5 (Jij) , 



PH m = ^5 {H, - Hr) + [Hi + Hr) , 



(2) 



(3) 



where J, c, and Hr are positive parameters. This distribution of exchange 
interactions is supposed to mimic a lattice of mean finite connectivity c. 

Using the replica method, it is not difficult to write the variational free 
energy [TT] 



r=Q («i,...,ar) 



In 



G(K})+/5i7^(T, 



(4) 



dHpH (H) Tifj exp 

where (3 = 1/ (ksT), T is the temperature, n is the number of replicas, 

br = cosh" {(3 J) tanh^ {(3 J) , (5) 

and the trace is taken over the set of replica spin variables {o"o}. The pa- 
rameter qai,...,ar is the expected value of the product aai---crar- The global 
order parameter G ({cTq}) is defined as 



G{{(ra})=C^ ^ brQo 



(6) 



In the n ^ limit, the minimization of this variational free energy with 
respect to the set of variables {qai,...,ar} leads to the stationary conditions 



Qai,...,ar 

where 



1 



+00 



dHpH (H) Tr^(Ja,...cr„^ exp 



G{{a^})+(3Hj2^c 



a=l 



dHpH (H) Tr^. exp 



a=l 



(7) 



(8) 



In the context of the replica-symmetric Ansatz, we define an effective field 
h, associated with an effective probability distribution p (h) , which is equally 
applied on all of the replica spin variables. We then write 



dai,...,ar 



dhp (h) tanV ipj) 



from which we have 



p{h) 



dHpH {H) 



dy_ 

277 



exp 



-iy {h-H) + G 



(9) 



• (10) 



Depending on the quantity of interest, it may be more convenient to work 
either with the global order parameter G or with the effective probability 
distribution p. In the replica-symmetric context, it is easy to see that 



G{{(y^})=G[Y,^A. 



(11) 



Therefore, taking into account the double-delta distribution Q, the extrem- 
ization of the variational free energy is reduced to the problem of searching 
the solutions of the non-linear integral equation 



p{h) 



+00 



dy_ 

271 



exp 



-iyh + In cosh {iyHn) — c + G 



[12) 



where 



/+00 
dxp (x) exp |y tanh^^ [tanh {f3J) tanh {l3x)] } . (13) 
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In the context of the rephca-symmetric Ansatz, it is known that the 
Viana-Bray spin-glass problem, with a symmetric distribution of exchange 
interactions and no external fields, can also be formulated in terms of a 
similar integral equation for the distribution of the effective fields jH] JO] ■ la- 
this spin-glass case, if we restrict to the analysis of the ground state (/9 — *• oc), 
it is possible to write an analytic solution as a sum of delta functions peaked 
at integer multiples of the variance J of the distribution of exchanges. In the 
present case, however, a similar solution requires the additional assumption 
that Hji/ J is restricted to the set of integer numbers. Under these conditions, 
in the ground state, the integral equation (fT^ can be exactly solved in terms 
of sums of delta functions. 

In the ground state (/3 oo), it is possible to show that 



+J 



dxp (x) exp {—iyj) + c dxp (x) exp (ixy) + 



+c / dxp (x) exp (iyJ) . 
J+J 

Assuming that Hji = uJ, with u; = 0, 1, 2, we can also write 



(14) 



G 



A + B exp (iyJ) + C exp {—iyj) , 



(15) 



with 



c = A + B + C, 



(16) 



and 



1a = ie(^-^) 
c 2 



B 



B 



L {2^/BC 



(17) 



C{l-t) 
B 



(2^BC{l-t)) + 



+ 



c{i-t) 

B 



' W {2^BC{l-t))}, 



(18) 
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for uj > 1, where I^^ (x) is the modified Bessel function. Both ferromagnetic 
{B 7^ C) and paramagnetic {B = C) solutions are represented by this ex- 
pression. The effective probabihty distribution p is written as a sum of delta 
functions, 



where 



afc = - exp {A - c) 



P 



(h) = J2 {h - kJ) 



(19) 



k=—oo 



k-\-uj 

' h+^(2VBC 



2VBC 



(20) 



3 Analysis of stability 

The analysis of stability of the replica-symmetric solutions is based on the 
investigation of the eigenvalues of the Hessian matrix associated with the 
variational free energy, given by Eq. which can be rewritten as 



PflG] 



)-2n 



2 ^ ^ 2cbr 

r=0 (ai,...,ar) 



In Tr^r exp 



G ({a„}) + In cosh ( (3Hr ^ (Xo 



a=l 



We then write 



Tr, 



S'Pf [G] 



5G(K})5G({rJ) 
which can be cast in the form 



ViM) = A^({(T„}) 



(21) 



(22) 



(f ({sq}) = cATr,- exp 



a=l 



+—TTr exp 



G (?) + In cosh {(3Hrt) + (^J^ t^Sc 
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-^G (s) Tir exp [G (?) + In cosh (^Hr?)] ^ ({rj) 



(23) 



where 

n n 

^=X1^"' s = ^s„. (24) 

a=l a=l 

In the n — i> hmit, it is easy to show that there is a constant eigenvector, 
if ({sq,}) = constant, with 1/c as the associated eigenvalue. 

According to the work of De Dominicis and Mottishaw ^T] , in the con- 
text of the rephca-symmetric approximation the space of 2" eigenvectors is 
spanned by a set of eigenvectors parameterized by functions of two variables, 
of the form 

V = <^{f,^} (S=; qa^,) , (25) 

where 



a=l 



a=l 



(26) 



and the spin configuration {/ia} is used to label the eigenvectors. From Eq. 
(|23|) ■ in the n ^ limit, we derive an integral equation for the eigenvalues. 



(3^' y) 



Gix] 



exp 

+ 00 



exp [c) 



im\ 



dmdr 



dudv 



G ( —J — imu — irv + In cosh (iHjim 



cosh {Pu + j3v) 
cosh {[3u — f3v) 

( 



m/2 



im ir 



+ 



dmdr 



dudv 



> (27r) 
X exp {imu + irv 



X + exp 



— J + In cosh {iHjim) — c 



cosh {(3J + (3u + (3v) 
cosh {pj — Pu — Pv) 

X cosh {PJ + Pu + Pv) cosh {pj - Pu - Pv)] ^ — — ^ ^ 

[cosh [p.J — Pu + Pv) 

X [cosh {pj + Pu- pv) cosh {pj - pu + pv)]~~ 



x — y 
4 



(27) 
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First, we find the longitudinal eigenvalues, in other words, the eigenvalues 
in the subspace spanned by (S"; qo-^) = ip (d). In the ground state, it is 
not difficult to see that these longitudinal eigenvectors can be written as 

<^L -T = Al + Bl exp {ixJ) + Cl exp {-ixJ) . (28) 



The problem reduces to the calculation of the eigenvalues of a 3 x 3 matrix, 
which are given by Ai = 1/c, associated with the constant eigenvector, and 

A2,3 = -il-A)±D, (29) 
c 

where A = cao is given by Eq. (fTTj). and D = (aia_i)^''^ is given by the 
expression 

D = lexpiA-c){ (^^^ J^+i (2v^) (2v^) + /^^ (2VBC 

with the coefficients A, B, and C, given by Eqs. (fT^ -()18 p . These longitudi- 
nal eigenvalues, however, lead to familiar mean-field results. At small values 
of the random field, the paramagnetic solution is unstable, while a ferromag- 
netic solution is stable. At large values of the random field, there is only 
a (stable) paramagnetic solution as in the case of spin glasses on the Bethe 
lattice (note that, at sufficiently large random fields, that is, for Hji — *■ 00, 
we have A,D —>■ 0). The critical border separating these paramagnetic and 
ferromagnetic phase is of order cJ (see Table 1 for specific values ). 

We now turn to the eigenvalues associated with the transversal sector. 
For /i = 0, and an eigenvector of the form 

(Pfi=o -75-7 = ^0 + Bi+ exp (ixJ) + exp {-ixJ) + 



(30) 



+B2+ exp (iyJ) + B2- exp {—iyj) + C++ exp {ixJ + iyJ) + 



+C+_ exp {ixJ — iyJ) + C_+ exp {—ixJ + iyJ) + C exp {—ixJ — iyJ) , 

(31) 
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the problem is reduced to the calculation of the eigenvalues of a 9 x 9 ma- 
trix. The nine eigenvalues of this transversal sector are given by Xti = l/c, 
associated with a constant eigenvector, 

At2,T3 = -{l-A), (32) 



At4,T5,T6 = -{1-A) + D, (33) 

and 

At7,T8,T9 = -{1-A)-D, (34) 

c 

which should be compared with Eq. for the non-trivial eigenvalues of 
the longitudinal sector. In contrast to the spin-glass case, the eigenvalues of 
this transverse sector do not lead to any additional instability. According 
to a numerical analysis of these eigenvalues, the replica-symmetric param- 
agnetic solution remains stable for sufficiently large random fields and the 
ferromagnetic solution is stable in its region of existence (see Table 1). 

The analysis of the transversal sector with fj. ^ can be carried out with 
same Ansatz, 

(f ' f ) = E {I'^^j) ^^--0 I) . (35) 

Although the secular matrix becomes infinite, it is easy to see that the eigen- 
values are still given by the same expressions of Eqs. (|32 p -(|34 p . Again, 
we conclude that the replica-symmetric solution is stable in the presence of 
sufficiently large random fields. 

In Table 1, we list the numerical solutions for the smallest eigenvalue, 
given by equation (jHH), with c = 3,4,5, and 6, for the paramagnetic (pm) 
and ferromagnetic (fm) solutions. There is no simultaneous instability of 
both solutions and thus no indication of breaking of replica symmetry. Note 
that the dashes in this table correspond to the absence of a ferromagnetic 
solution (in which case the paramagnetic solution is stable). 



Table 1 
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Hr/J 







1 




2 


3 


4 


c 


Xpm 




-0.1798 




■0.0661 


0.1030 


0.2253 


0.2919 


3 


Xfm 


0, 


.2738 


0. 


,1157 








3 


Xpm 




■0.1876 




■0.1159 


0.0129 


0.1212 


0.1906 


4 


Xfm 


0, 


.2302 


0, 


.1669 








4 


Xpm 




-0.1876 




■0.1377 


-0.0369 


0.0578 


0.1258 


5 


Xfm 


0, 


.1930 


0, 


.1666 


0.0819 






5 


Xpm 




■0.1847 




■0.1476 


-0.0667 


0.0161 


0.0808 


6 


Xfm 


0, 


.1642 


0, 


.1529 


0.1153 
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4 Conclusions 

We have investigated the stabihty of the rephca-symmetric solutions of a 
random-field Ising ferromagnet on a lattice of finite mean connectivity. At 
low temperatures and for sufficiently large random fields, the analysis of the 
eigenvalues of the Hessian matrix associated with the variational free energy 
leads to stable replica-symmetric paramagnetic solutions (at smaller random 
fields, the replica-symmetric ferromagnetic solution is stable). The present 
calculations do not support the existence of a glassy phase, as suggested by 
earlier proposals However, a more detailed analysis of the phase 

diagram, in terms of field and temperature, still demands considerable work, 
including both analytical and refined numerical calculations. The assumption 
of a discrete distribution of effective fields, which was written as a sum of 
delta functions, may not capture the subtleties of the glassy behavior. As in 
the spin-glass case, we cannot rule out the existence of field-induced glassy 
and mixed ferromagnetic-glassy phases. 

It is interesting to point out a connection with the recent work by Pastor 
and collaborators A truncation of the variational free energy, given by 
Eq. (jD), leads to a "high-temperature approximation," 

1 X c c 
fapp = 15 -{o ^^^^ ^l^J) yi^l+n [^^^^ *^^'^)]^ 5Z ^"/3~ 

' a a</3 
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-In 



app 



(36) 



where 



Zapp 

Tr^ exp{ctanh (/5J) ^ cra+ 



+c [tanh ipJ)f ^ qaf3(^a(^p + pH^^ a^}. (37) 

a</3 a=l 

The results of Pastor et al. 3J for the paramagnetic phase are recovered if 
we introduce the additional approximation tanh/? J = PJ + O [(/3J)^], and 
discard higher-order terms. In this approximation, for a Gaussian distribu- 
tion of random fields with variance Hji, the replica-symmetric paramagnetic 
phase is unstable along the field axis. 
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